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ABSTRACT 

Massive galaxies with a half-mass radius 1 kpc are observed in the early universe 
[z 2), but not in the local universe. In the local universe similar- mass (within a 
factor of two) galaxies tend to be a factor of 4 to 5 larger. Dry minor mergers are 
known to drive the evolution of the size of a galaxy without much increasing the 
mass, but it is unclear if the growth in size is sufficient to explain the observations. 
We test the hypothesis that galaxies grow through dry minor mergers by simulating 
merging galaxies with mass ratios of q =1:1 (equal mass) to q =1:160. In our N- 
body simulations the total mass of the parent galaxy doubles. We confirm that major 
mergers do not cause a sufficient growth in size. The observation can be explained 
with mergers with a mass ratio of q =1:5-1:10. Smaller mass ratios cause a more 
dramatic growth in size, up to a factor of ~ 17 for mergers with a mass ratio of 1:80. 
For relatively massive minor mergers q ^ 1:20 the mass of the incoming child galaxies 
tend to settle in the halo of the parent galaxy. This is caused by the tidal stripping of 
the child galaxies by the time they enter the central portion of the parent. When the 
accretion of minor galaxies becomes more continuous, when q ^ 1:40, the foreign mass 
tends to concentrate more in the central region of the parent galaxy. We speculate that 
this is caused by dynamic interactions between the child galaxies inside the merger 
remnant and the longer merging times when the difference in mass is larger. These 
interactions cause dynamical heating which results in accretion of mass inside the 
galaxy core and a reduction of the parent's circular velocity and density. 

Key words: galaxies: evolution - galaxies: interactions - galaxies: kinematics and 
dynamics - methods: numerical 



1 INTRODUCTION 



Hierarchical structure formation drives t he growth of 
mass and size of galaxies with redshift (Peebled Il98d : 
lLacev fc Cold [l993f ) . In this picture the mass and size of 
galaxies grow at a similar rate. Recently a population of 
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mass, which cannot be explained from major mergers as 
obse rved in hierarchic al structure f ormat ion simulations 
(e.g. iNaab et all (|200rj '): lOser et all (|2012r i and references 
therein). 

It has been suggested that dry minor mergers can cause 
a considerable growth in size without much increasing the 
mass of the galaxy, whereas major mergers tend to increase 
the mass without much increasing the size IIMiller fc S mith 
198d: INaab et alj|2005 lOser et al.|[201ol . |2012l ; ITrujillo et all 



However at low redshift t hese g alaxies are much 
more 
20091: 



20111 ). 
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dTruiillo et al l 12009] : 



Valentinuzzi et all l20ld : ITrujillo et all I2012T ). whereas 



there is a rich population of elliptical galaxies with 
similar mass but which are considerably larger in size 
l|Martinez-Manso et al.l [2011). This suggests that small 
and massive galaxies grow in size without acquiring much 
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Several studies on the topic seem to contradict each 
other, some argue that minor mergers can drive the ob- 
served growth e.g. dBezanson et al.ll2009l : INaab et al. I l2009l : 
Hopk ms et al.ll2009l . l20ld 'l. Others claim that the observed 
size growth in simulations is not suffici ent if one takes into 
account cosmological scal ing relations (|Nipoti et al. 2009a, 
120121 : ICimatti et al.112012] ). Encouraged by this discrepancy 
in the literature we decided to perform a series of simula- 
tions in which we model the encounters between compact 
massive galaxies and smaller, lower-mass counterparts. In 
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our parameter search we included mass ratios from 1:1 all 
the way to 1:160. Even t hough such low-m ass encounters are 
unlikely to be common l|Oser et alJl2012h . they approach a 
continuous in-fall of material on the parent galaxy. 

With the merger simulations we study the growth in 
mass, size and the effect on the shape of the merger prod- 
uct. Our simulations are performed using the Bonsai tree 
code (|Bedorf et al . 2012) with up to 17.2 million equal- 
mass particles. This includes the dark-matter as well as 
the baryonic matter. In our simulations we recognize two 
regimes of galaxy growth. For q ^ 1:20 we confi rm the inside- 
out growth, as discussed by ()Hilz et al.ll2012bh . whereas for 
q 1:40 galaxies tend to form outside-in, meaning that the 
mass is accreted in the core of the primary galaxy rather 
than accreted on the outside. This change in behavior is 
caused by the self-interactions of minor galaxies for mass 
ratios ,$ 1 : 40. 

Recently jHilz et all 1201281 ) studied the size increase 
from major mergers down to 1:10 minor mergers. Their re- 
sults are consistent with our findings for compar able simu- 
lation parameters. The major difference between |Hilz et al.l 
l2012al ) and the results presented here are the details regard- 
ing the orbital parameters of the merging galaxies; they drop 
their minor galaxies in one-by-one, whereas we initialize all 
galaxies in a spherical distribution at the start of the sim- 
ulation. As a consequence in some of our simulations the 
merging process with many minor mergers is not completed 
by T = 10 Gyr. While the results of our 1:1 to 1:10 mergers 
ar e consisten t with t he gro wth observed in the s imulations 
of iHilz et all j2012allbT) and lOogi fc Habd |2012l ). they are 
inconsistent with Nipot i et al. I (|2009bl ). The reason for this 
discrepancy cannot be the st eeper dens i ty slope of the stellar 
bulge, as was suggested by IHilz et al.l |2012al ). because we 
tried even higher density slopes for the minor child galax- 
ies and found that this does not affect the growth of the 
merger product (see Ap pendixlBl). The discre pancy between 
our results and those of lNipoti et all l|2009bl laT) could be ex- 
plained by the a veraging of the sim ulation results, as was 
also suggested bv lHilz et al.l (|2012al ). 



2 CONSTRAINING THE MODEL 
PARAMETERS 

To perform the simulations we use an updat ed version of 
the g ravitational TV-body tree-code Bonsai (|Bedorf et al.l 
l2012h . All simulations are performed on workstations with 
NVIDIA GTX480 graphical processing units. Simulations 
are run with 1 to 18 GPUs, depending on N. 

The code has three important parameters; the choice 
of the time step dt, the softening length e and the tree- 
code opening angle 9. The opening angle is set to 9 = 0.5 
whereby we use the min imum distance opening criterion 
llSalmon fc Warren! QUI ), which is computation ally more 
expen sive than the improved Barnes-Hut method (|Dubinskil 
1996) by about a f actor of two, but offers better accuracy 
i|Bedorf et al]|2012l ). 

After having fixed 9, the choice of dt and e are quite 
critical for the quality of our results. We measure the qual- 
ity of the simulation based on the error in the energy and 
the change in size of a galaxy model in isolation. We de- 
termine the optimal values for dt and e by performing an 



analysis in which we run a series of models using a range 
of values for dt and e. Although the time step and soft- 
ening are discussed extensively in prev i ous literature (e.g. 
|Merrittlll996l : lAthanassoula et aLll2000l : lDehnenll200ll )) we 
chose to do extra tests our-self. This for a couple of rea- 
sons, first most of the previous discussed methods propose 
a softening that scales with the number of particles. How- 
ever, we use different number of particles for the primary 
galaxy and child galaxies so we had to find a configuration 
that works for both. Second, our tree-code uses an accurate 
multipole expansion, in combination with a opening angle 
criteria that is different from used in the previously men- 
tioned papers. And finally we use two component models 
with power law density distributions while most examples in 
the previous mentioned papers are based on one component 
spherical Plummer models. Therefore we decided to rather 
empirically test how the softening and time-step affected the 
properties of our initial conditions instead of taking over a 
previously result obtained in a different setting. The initial 
conditions for these isolated galaxies are generated using 
GalactlCS ilKuiiken fc Dubinskilll995l : IWidrow fc Dubinskil 
120051 : IWidrow et al.ll2008l ) andconsist of a dark matter halo 
and a baryonic bulge. All models represent elliptical E0 
galaxies. The bulges are modeled with a Sersic index of 3 
and the dark-matter halo is represented by an NFW profile 
l|Navarro et al.ll 19961 ). 

In our simulations we collide the primary galaxy with a 
number of child galaxies, which are smaller and less massive. 
The mass ratio and the number of child galaxies which inter- 
act with the primary is one of the free parameters in our sim- 
ulations. We generate the child galaxies using GalactlCS by 
adopting the appropriate mass and size of the dark-matter 
halo both of which we assume to be a constant factor smaller 
than in the primary galaxy. This factor is equal to the mass 
ratio of the child, e.g. for a 1:10 child the mass and cut-off 
radius is 10 times lower than that of the primary galaxy. 
As a consequence the density of the child galaxy is smaller 
by a factor ~10 than that of the primary. In AppendixlBl 
we demonstrate that the effect of our selected size and mass 
(and consequently density) of the child galaxy is not critical 
for our conclusions. In Fig. [T] we present the density profiles 
of the baryonic bulge and dark matter halo of our selected 
initial conditions at zero-age, at T = 1 Gyr which is the mo- 
ment we intend to use them for the merger simulations in 
§13] and at T — 10 Gyr. It takes a few crossing times for the 
galaxies to relax after generation, therefore we first let them 
evolve for T = 1 Gyr before we use them in a merger config- 
uration. The differences in the density profile for the same 
component are hardly noticeable as a function of time. To 
prevent spurious mass segregation all particles in one simula- 
tion have the same mass. Mass differences between galaxies 
and baryonic and dark matter particle distributions are cre- 
ated by using different amounts of particles for each galaxy 
and galaxy components (see Tab. [T] for details). 

The isolated galaxy simulations were performed for a 
galaxy with a total mass M parcnt = 2.2 x 10 12 Mq, and a 
half-mass radius of 1.36 kpc. which we call the parent or pri- 
mary galaxy and for each of the child galaxies. The child 
galaxies have a mass that is q times lower than that of 
the primary. For example the q =1:10 child has a mass of 
M chi id = 2.2 x 10 11 M Q . To make sure that the child galax- 
ies consist of enough particles to be accurately resolved and 
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Figure 1. Density profiles of the HR baryonic bulge and dark- 
matter halos for the primary and child galaxies, which have 
evolved in isolation. For each component we present the density 
profile at zero age, at an age of 1 Gyr and at 10 Gyr. However, be- 
cause the lines are almost indistinguishable and to prevent clutter 
we show each epoch with the same line style. The expansion in 
time of the child halo is caused by the escape of a few dark-matter 
particles. This is also the main reason for the expansion of the 
selected model, noticeable in the right hand-panel of Fig. [3] The 
vertical line near R = 0.05 kpc indicates the softening length in 
the HR, simulation. 



stay stable in isolation we use different resolution models 
for the initial conditions. Depending on the number of child 
galaxies (mass ratio) used. We found that using a minimum 
number of 10 3 (10 4 ) particles for the bulge (dark matter) 
component proofed to be enough to keep the galaxies stable 
(Tab. [[]) . 

We refer to low-resolution simulations (LR) for those in 
which the primary galaxy was modeled with TV = 2.2 x 10 5 
and the 1:10 child with TV = 2.2 x 10 4 . For the high- 
resolution (HR) simulations we adopted TV = 2.2 x 10 6 and 
TV = 2.2 x 10 5 for the primary and 1:10 child galaxies, re- 
spectively. The child galaxies with a different mass ratio are 
scaled in a similar way as the 1:10 child, namely the num- 
ber of particles of the primary divided by the mass ratio. 
Each model was run for 10 Gyr. In Figs.[2]and[3]we present 
the results of the HR simulations; the error in the energy 
AE = (Eiooyr — Eq)/Eo for the left-hand panels and the 
expansion factor of the bulge, 7(t), in isolation in the right- 
hand panels. 



7(t) = Rh-t/Rh-0Gyr 



(1) 



Here Ru-t is the half-mass radius of the baryonic component 
of the galaxy at time t. In Fig. [2] we present the result for 
the primaries, and Fig. [3] for the 1:10 children. 

As expected, AE decreases for decreasing dt and in- 
creasing e, for the primary as well as for the child galaxy. 
The growth in size of the galaxies becomes smaller with de- 
creasing dt. However, for e this is not a trivial matter. For 
too large a softening the galaxy tends to expand substan- 
tially, but the effect is not as dramatic as for a too small 
value of e. The growth for large e is caused by spurious sup- 
pression of relaxation in the large e runs. The growth for 
small e is numerical, caused by strong encounters which are 
not resolved accurately in our simulation. This limitation in 



our numerical solver is also noticeable in the energy error 
for small e (see Figs. [2] and Fig. [3]). 

Because we intend to study the expansion of the pri- 
mary galaxies caused by mergers they should not expand in 
isolation. For our production simulations (see §13]) we select 
the e which show the least expansion when the galaxy is run 
in isolation. The time step is chosen as large as possible (for 
performance) but giving the smallest possible energy error. 
We make the choice of dt and e based on the simulations 
for the primary galaxy, which for the child also turns out to 
be the bests choice. For the HR simulations dt — 1/192 and 
e = 0.05, and for the LR simulations dt — 1/128 and e = 0.1. 
Using these settings the energy error of all the merger sim- 
ulations stays below <C 10 -4 . These p arameters are pre- 
sente d in dimension-less TV-body units l|Heggie fc Mathieul 
1986), but in physical units they translate to a time step 
of dt ~ 0.078 Myr for the LR and dt ~ 0.052 Myr for the 
HR simulations, and a softening length of e ~ 100 pc, and 
e ~ 50 pc for the LR and HR runs, respectively. 

Comparing the most optimal settings with the re- 
sults of previous work on softening va lues (|Merrittj 1 19961 : 
lAthanassoula et al.l 2000; De hnenl l200lh . The chosen soft- 
ening values are compa rable to the valu es advised in 
l|Athanassoula et alj|200ch and (lDehnenll200 ifi with our cho- 
sen values being slightly larger to keep the energy error less 
than 10~ 3 while the size growth sta ys smalle r than 10%. The 
difference between our results and iMerrittJ (|1996l ) is some- 
what larger, but still within a factor of three. The difference 
can be attributed to the fact that we use more particles and 
use a tree-code instead o f direc t summation (see also Table 
1. of lAthanassoula et all (Il998l )). 

Most of the expansion of our galaxies occurs in the first 
few 100 Myr of the simulation, and therefore all galaxies are 
evolved to an age of 1 Gyr in isolation before they are used 
as initial galaxies in the merger simulations. 

As as extra test we ran several of t he isolated galaxies 
using the direct TV-body code phi GRAPE (jHarfst et al.ll2007h 
to confirm that the observed AE and expansion are not the 
result of using an approximation based simulation code. The 
results of these simulations are consistent with the results 
of Bonsai, but show a AE of ~ 10~ 6 , because of the higher 
accuracy of phi GRAPE . 

With this choice of initial conditions, the child galaxies 
in the low resolution simulations still expand by at most a 
factor of 2 before the merger starts. To check how this affects 
our results we redo the simulation but at the start the child 
galaxies are represented as point masses. The point mass 
is replaced by a child galaxy as soon as it enters the dark 
matter halo of the primary galaxy. This simulation gives a 
nearly identical expansion of the primary galaxy to that in 
which the child was resolved from the start of the simulation. 



3 INITIALIZING THE GALAXY MERGERS 

We use the isolated galaxies discussed in the previous § to 
study the effect of mergers. For this purpose we adopt the 
primary galaxy and have it interact with an identical copy or 
with a number of child galaxies. We refer to minor mergers 
when the primaryxhild mass ^ 1:5. Major mergers in our 
study have equal mass. 
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Figure 2. Results of the test simulations for the HR primary galaxies. The left panel gives the energy error AE, the right panel gives 
the size increase 7 after lOGyr (see Eq.(2j- Both panels present the information as a function of the time step (dt) and the softening 
length (e). The simulations are performed at the grid-points in dt and e, as indicated along the axes. The color scheme shows a linear 
interpolation between the grid points. The scaling to the colors are provided along the right-hand side of each panel, and exhibits a 
log-scale for the left hand panel and a linear scaling for the right-hand panel. 
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Figure 3. Results of the test simulations for the HR child galaxies, which are ten times less massive than the primary. See Fig. [2] for a 
further description of the figure. 



3.1 Configuring the major mergers 

In the major merger simulations two identical copies of the 
primary galaxy are placed on an elliptical orbit. The initial 
distance between the two galaxies is 400 kpc, which exceeds 
the size of the dark-matter halos, which is about 200 kpc. 
The relative velocity is chosen such that the minimal dis- 



tance of approach during the first perigalactic passage P 
varies. We vary P from head-on (Okpc) to 10, 50, 75 and 
100 kpc for the widest encounters. All our merger simula- 
tions are run until 10 Gyr, after which both galaxies have 
merged to a single galaxy, except for the P = 100 kpc con- 
figurations (see § [4j . 

Because we use two identical copies of the primary 
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Table 1. Galaxy properties. Characteristics of the base models used in the simulations. The first column indicates if the galaxy is the 
primary galaxy or one of the minor merger child galaxies. The second column indicates if the properties are for the dark matter halo 
or for the baryonic bulge and directly after in the third column the mass of the galaxy component. The next three columns show the 
properties when using low resolution models, the number of particles, mass per particle and the softening used. The final three columns 
show the number of particles, mass per particle and softening used when using high resolution models. 



galaxy the galaxies will be co-rotating during the collision. 
To study the influence of the rotation angle on the size of 
the merger product we run the simulations for P = and 
P — 50kpc with co-rotating and counter- rotating galaxies. 
This is done by changing the inclination and Q of one of 
the two galaxies between ° and 270 ° in steps of 90 ° using 
a rotation matrix. After the rotation the galaxies are put 
on their elliptical orbit as described in the previous §. The 
stability of the isolated galaxies is not affected, because the 
positions and velocities are rotated in a consistent way. We 
did not perform this study on the other major mergers, be- 
cause the effect of the inclination and il on the expansion 
factor turns out to be negligible (variation of the final size 
between the 16 configurations is ~ 0.02 kpc). 



3.2 Configuring the minor mergers 

The easiest way to generate an initial configuration with 1 
primary galaxy and N child galaxies is by randomly posi- 
tioning them in a spherical distribution. For each simula - 
tion we generate a Plummer distribution l)PlummerJ 1 19151 ) 
of TV + 1 particles. We then select a random particle and 
reposition it in the center of mass. We assign the mass of 
the parent galaxy to this central particle. The other parti- 
cles are assigned the mass of the child galaxies. The total 
mass of all child galaxies in each simulation is the same as 
that of the primary galaxy. The model is then rescaled to 
a virial radius, R, and virial ratio (temperature), Q. We 
run 10 sets of initial conditions for each selected mass ratio. 
In the next section we discuss the results of a number of 
simulations with varying mass-ratio, R and Q. Using this 
procedure we created a range of minor merger models with 
varying R, Q — and N, see Tab.(2] In Section 13.2.11 we 
relax this assumption and allow Q > to study the effect of 
the temperature on the size of the merger product. 
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Table 2. Initial conditions for the minor merger simulations. The 
first column indicates the used mass ratio. The second and third 
column indicate the virial radius, R, and virial temperature, Q, 
to which the configuration is scaled. The fourth column indicates 
the used resolution, either low-resolution (N = 4.4 X 10 5 ) or high- 
resolution (N = 4.4 X 10 6 ). The final column indicates the name 
used in the text when we refer to the configuration. 

For the 1:10 minor mergers we adopted R — 200 kpc as 
well as R — 500 kpc. In the R = 200 configurations the child 
galaxies have a ~ 50% probability of being placed initially 
within the dark matter halo of the primary. For the R = 
500 configurations, which are based on the same random 



based on same Plummer distribution as the R = 200 models. 
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Plummer realization but scaled to a larger virial radius, this 
is only ~ 15%. 



3.2.1 Plummer models - Varying virial temperatures 

In order to study the effect of Q on our results we selected 
one of the 1:10 models. The merger remnant of the selected 
model has roughly the same mass in both the R = 200 and 
R — 500 configurations indicating that they were involved 
in an equal number of mergers by the end of the simulation. 
We used this initial configuration to vary Q between (cold) 
and 1 (warm), see Tab. [5] Each run was repeated with R = 
200 kpc and for R — 500 kpc. This results in a total of 12 
different configurations, all of which are based on the same 
spatial distribution but with different virial radii and virial 
temperatures. 



4 RESULTS 

During the simulation snapshots are taken every 100M year 
and post-processed using tipsy 0. We compute the density 
and cumulative mass profiles which are subsequently used 
to determine the half-mass radius (Rh), measured using the 
baryonic particles only. To compute the profiles tipsy com- 
putes the densest point in the simulation and then uses 
spherical shells centered on the densest point to bin the par- 
ticles. To compute Rh we can not simply take the radius that 
contains half of the total mass in the merger system, because 
depending on the time of the snapshot, not all child galaxies 
have merged with the primary galaxy. Therefore we analyze 
the cumulative mass profile (based on all baryonic particles 
in the system) of the galaxy and take as total mass the point 
where the profile flattens. This flattening indicates that all 
child galaxies within that radius have merge with the pri- 
mary galaxy and that the unmerged child galaxies are too 
far away to be counted as part of the merger remnant. This 
total mass is then used to determine Rh- We illustrate this 
procedure in Fig. [4] 

In Fig.[4]we present the cumulative mass-profiles of one 
of the 1:10 mergers with R = 500 and Q = 0, for the ma- 
jor merger with P = 50 kpc and the profile of the isolated 
primary galaxy. During the first four Gyr of evolution the 
cumulative mass of the minor merger simulation gradually 
increases each time one of the child galaxies is accreted. We 
determine the total mass of the merger remnant by deter- 
mining the first moment that the cumulative mass profile 
becomes flat. This happens, for example, for the 1 Gyr sit- 
uation at a mass of M ~ 2.2 x 10 11 Mq, and for 2 Gyr at 
M ~ 3.2 x 10 11 Mq. We further notice that the minor merger 
remnant has grown in size with respect to the major merger 
remnant, even though it has accreted less mass. The major 
merger remnant ends up with twice the mass of the primary 
galaxy, e.g. a complete merger of all mass in the system. 
Whereas in the minor merger only 8 child galaxies have been 
accreted in the 4 Gyr after the start of the simulation. At 



2 http:/ /www-hpcc. astro. washington.edu/tools/tipsy/tipsy.html 

3 We chose this method over a method where we check if particles 
are bound, since in this method we include child galaxies that 
would add to the luminosity of the remnant galaxy if they would 
be observed. 
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Figure 4. Cumulative mass profiles of the bulge. Shown are the 
profiles of one of the 1:10 merger simulations with R = 500 and 
Q = and the major merger simulation with P = 50 kpc. The 
thick solid line shows the cumulative mass profile of the main 
galaxy at the start of the simulation. The horizontal lines show 
the half-mass of the major merger (top line at 2.0Mq) and of 
the minor merger (bottom line at 1.86Mq). The remaining lines 
show cumulative mass profiles of the minor merger as indicated. 



the end of the simulation (T=10Gyr) 9 child galaxies have 
merged with the primary galaxy. The last unmerged child 
appears in the cumulative mass distribution (striped line 
with solid circles), represented by the bump in the line at 
r ~ 1800 kpc (Fig. |4}. As a consequence the total mass of 
the merger product at the end of the simulation is larger for 
the major merger than for the minor merger. 



4.1 The growth of the primary due to subsequent 
mergers 

We run each major merger model 16 times (with different in- 
clination and Q ) and average the Rh . For each set of minor 
merger initial conditions (R and N, see Tab.0) we perform 
10 simulations and as in the major mergers we average Rh- 
These results are presented in Fig. [5] for the growth in size, 
and Fig. [5] for the growth in mass. The shaded areas are 
the ler deviation from the mean, but for the major mergers 
the areas are barely visible in comparison with the width 
of the line. This indicates that the inclination and the an- 
gle fl hardly affects the growth of the merger product. Only 
during the merger event itself some deviation is noticeable 
(at about 3 Gyr for the simulation with P = 50 kpc) . The 
head-on configuration (P = 0) grows from an initial size of 
Rh ~ 1.3 kpc to 3.4 kpc, whereas the P = 10 to 75 kpc sim- 
ulations grow to only about 2.4 kpc (only the P = 50 kpc is 
presented in Fig. [5}. The P = 100 kpc did not experience a 
merger in our simulation within the time frame of 10 Gyr. 
For the minor merger configurations the result is quite dif- 
ferent as can be seen in Fig. [5] The shaded areas indicate 
that there is variation in size between the different configu- 
rations. The size growth however is larger than that of the 
major mergers. The results indicate that the growth in size 
is directly related to the mass ratio of the child galaxies, the 
larger the difference in mass the larger the growth in size. 
For the 1:10 simulations the Rh is at least 7 kpc and for the 
1:80 simulations the average Rh is at least 20 kpc while at 
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Figure 5. Evolution of Rh- The plot shows the mean (solid and 
dotted lines) and the standard deviation (shaded areas) of the Rh 
based on different merger configurations. The solid lines show the 
results of the minor merger configurations. The dotted lines show 
the major merger configurations with P = and P = 50. The 
minor mergers are each averaged over 10 different Plummer real- 
izations, the major mergers are averaged over 16 different combi- 
nations of Euler angles. For 1:80 we filtered out size jumps caused 
by invalid size detection. 
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Figure 6. Evolution of the merger products mass. The same 
simulations are shown as in Fig. [5] but instead of the radius the 
average mass of each set of simulations is shown. 



the start of the simulation the primary galaxy has a Rh of 
1.3 kpc. 

To study the effect of increasing the number of minor 
mergers (and decreasing the mass per child galaxy) in more 
detail we start by comparing the major merger with the 1:5 
minor merger. This case results in a Rh of 5.1 kpc, averaged 
over 10 simulations. The lcr deviation indicates that the re- 
sult is rather consistent when rerunning with a different ini- 
tial Plummer realization. The extreme values of Rh after 
10 Gyr range between 4 and 6 kpc. Even the increase due to 
the 1:5 minor merger exceeds the P = major merger case 
by a factor of ~ 1.2 to 1.8. By the end of the simulation at 
10 Gyr on average 4 out of the 5 minor galaxies have merged 
with the primary galaxy, and we expect that if the last child 
galaxy would be accreted the increase in size would be even 
larger. The 1:10, R200kpc minor mergers show an average 



size of 7.6 kpc, and 8.5 kpc for the 1:10, R500kpc models; an 
increase of 5.8 to 6.5. The R = 200 kpc simulations merge 
on a shorter time scale compared to the R = 500 kpc, which 
results in a smaller expansion, but higher mass of the sim- 
ulations performed in a smaller volume. This is also visible 
in Fig. [H] where the R = 200 kpc increases faster in mass 
and ends up being more massive than the R = 500 kpc case. 
Even though R — 500 kpc results in a less massive merger 
remnant than the R — 200 kpc simulations, it grows to a 
larger size. This is caused by the disturbing effect that the 
dark-matter halo of the primary galaxy has on the child 
galaxies. With more galaxies placed outside the halo for the 
R — 500 kpc case this effect will be more pronounced. While 
we continue to increase the number of child galaxies the size 
of the merger remnant continues to grow. We stopped with 
N = 80 children, at which moment the increase in size com- 
pared to the major merger exceeded a factor 17. The main 
reason for the growth of the merger remnant when increasing 
the number of child galaxies is by the heating of the merger 
remnant. For each doubling of the number of child galaxies 
from 1:5 to 1:10, etc., the growth in size of the merger prod- 
uct is about constant by ~ 5 kpc. The growth in mass for 
each of these initial configurations is comparable, except for 
the 1:80, which grows less quickly. The smaller mass growth 
in the 1:80 is attributed to increased dynamical interactions 
among the child galaxies. These interactions cause the merg- 
ing times per galaxy to take longer than when they would 
merge sequentially. 

The combined growth in size and mass indicates an 
other effect, namely the larger the mass ratio the higher the 
growth in size. This is indicated in Fig. [7] where we present 
Rh as a function of the merger remnant mass (which is a 
function of time). This figure quantifies the growth in size 
of the merger remnant while increasing the number of merg- 
ers. Note that the mass ratio between the merger remnant 
and the child galaxies changes over time. For example the 
1:10 minor mergers start with a mass difference of a factor 
10 between one child galaxy and the primary galaxy, but 
after the primary galaxy has merged with 5 child galaxies 
the mass ratio between the merger remnant and one child 
galaxy is 1:15. Interestingly it is not necessarily the mass 
that drives the merger remnant size, but it is the number 
of children that merges with the primary galaxy. The solid 
curves are fits through the data, and represent exponen- 
tial curves. The dashed line is the observe d scaling relation 
of R h oc M 2 04 (|van Dokkum et al.|[201ol ). As we can see 
our minor mergers are above this relation while the major 
mergers are clearly smaller than the observational results. 
Indicating that the minor mergers can cause a size increase 
comparable or even larger than that observed. 

4.2 The effect on the galaxy shape due to 
subsequent mergers 

Now that we have established that the minor mergers can 
have a pronounced effect on the size of the galaxy remnant 
it will be interesting to see if we can recognize this galaxy 
growth by the shape of the merger remnant. To test the 
effect that the mergers have on the shape we measure the 
semi-principal axes a, b and c of the ellipsoidal with a >= 
b >= c. The axes lengths are obtained by computing the 
eigenvalues of the inertial tensor of the galaxy. Using this 
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Figure 7. The mass of the merger remnant, M, versus Rh for 
the different configurations. The solid lines are exponential fits 
through the results of the minor merger configurations. The major 
mergers are not fitted. The thick dashed li ne (bottom line) is 
the o bserved scaling relation of R h oc M 204 Jvan Dokku m et al.l 
l20ld) . Configurations are plotted after each merger event (but 
filtered to not show fly-by events) and therefore occur multiple 
times. The exponents of the fits are (from top to bottom), 0.9, 
0.84, 0.74, 0.62, 0.57 and 0.48. If we make an exponential fit to 
the observed scaling relation than this would have an exponent 
of 0.5. 



we plot the three axis ratios b/a, c/a and b/c, if these three 
ratios are 1 then the elliptical is perfectly spherical. If one 
of the axis is much smaller than the other two the galaxy is 
a flattened spherical. 

In Figs.[8]and [9]we present how the shape of the merger 
remnant (measured at Rh) evolves over time (baryonic com- 
ponent only). For the major mergers we used the default con- 
figurations with P — and 50 kpc and for the minor merger 
configurations we randomly took one of the 10 realisations 
for each of the used mass ratios. The x-axis shows the simu- 
lation time while each of the three panels presents one of the 
aforementioned axis ratios. Visible is that all galaxies start 
of as a perfect sphere, which is the initial condition, but over 
time the merger remnants deform. The major mergers show 
the highest degree of deformation with large differences be- 
tween the different axis. The effect is most severe for the 
heads-on collision which is expected as it is the most vio- 
lent collision of all our configurations. Although the major 
merger remnants become slightly more spherical after the 
mergers are complete the galaxies will not become perfectly 
spherical within a Hubble time. 

For 1:5 we notice, as with the the major mergers, a clear 
deformation from spherical, but the effect is less strong than 
for the major mergers. The individual merger events of 1:5 
are visible by the spikes in the axis ratios. The impact that 
a merger event has on the shape of the merger remnant 
becomes smaller when the mass difference between merger 
remnant and the infalling galaxy increases. For the 1:10 con- 
figurations, we do see slight deformations during the merger 
events (the shark tooth's in the lines), but when the merger 
is complete the remnant settles back to spherical. The same 
effect is visible for 1:20, 1:40 and 1:80 (Fig. [9]). The final 
merger remnant of the minor merger configurations is no 
longer perfectly spherical, but it clearly is not as flat as the 
major merger configurations. Indicating that the accreted 



mass is distributed differently through out the merger rem- 
nant in the minor merger simulations than in the major 
merger simulations. We cover the location of the accreted 
mass in more detail in Section [5] 

In Fig. [10] we present the galaxy shape at T = lOGyr 
from the core of the merger remnant up to 2 x Rh, while 
in the previous paragraph we looked at the shape of the 
merger over time and only at the Rh location. Because Rh 
is different for each configuration the lines representing the 
configurations in Fig. [10] have different lengths. The minor 
merger results are averaged over the different realisations of 
the model. We notice that the larger the differences in mass 
between the the primary galaxy and the child galaxies the 
more spherical the merger remnant stays (axis ratios near 
1). The differences between the major and minor mergers 
are especially visible on the outsides of the merger remnant 
(r > Rh)- The major mergers are much more deformed from 
perfectly spherical into flattened ellipticals while the minor 
mergers stay near spherical all the way to the outskirts of 
the galaxy. The difference in the core of 1:5 and 1:10 with 
respect to the other configurations is the result of running 
these simulations with only N = 4.4 x 10 5 particles and 
accompanying larger softening. We checked this by compar- 
ing the individual results of one of the major mergers and a 
1:10 high resolution simulation with the comparable config- 
urations in low resolution. In the high resolution simulations 
we do not observe the deformation in the core, instead the 
axis-ratios have a similar profile as the other high resolution 
configurations. 



4.3 The effect of the virial temperature 

In the previous section the virial temperature (Q) of the 
minor merger configurations was always set to zero. This 
resulted in cold collapse and consequently the fastest possi- 
ble merger. In this section we investigate the configurations 
described in Section [3.2.11 where Q is varied. The results are 
presented in Fig. [TTJ where we present Rh as a function of the 
merger remnant masfl (which is a function of time) . In the 
top panel the results for the R = 200 kpc are presented and 
the bottom panel presents the results for the R = 500 kpc 
configurations. The results for the major mergers with P — 
0, 10 and 50 kpc are presented in both panels. 

The major mergers are the same as the ones presented 
in the previous section and end up with a final mass that 
is double that of the isolated galaxy (M = ±4 x 10 11 Mq). 
Indicating that in all major merger configurations the two 
galaxies fully merge. The horizontal distribution in mass is 
caused by accretion of stars that were flung out and later fell 
back onto the merger product, thereby increasing the mass 
of the merger remnant. 

The merger events that occur in the minor merger sim- 
ulations can be seen in the horizontal distribution of the 
points. The horizontal distance between the points compares 
to a separation equal to the mass of a 1:10 child galaxy. 
Even though the minor merger configurations are based on 
the same spatial distribution, they are not per se involved 
in the same amount of mergers, this depends on R and Q. 



4 We filtered the results to not include fly-by events which tem- 
porarily increases Rh- 
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Figure 8. Axis ratios of the merger remnant. The three panels present the axis ratios b/a, c/a and b/c of the merger remnant (top to 
bottom, with a ^ b c). As location to measure the shape is taken the mean Rh of the configuration. If all axis ratios are equal to 1 
the shape is perfectly spherical. Shown are the results for the head-on major merger (P = Okpc), for the major merger with P = 50kpc, 
for a random 1:5 merger and for a random 1:10 merger once with R = 200 kpc and once with R = 500 kpc. 
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Figure 9. As Fig. [8] but now for one of the 1:20, 1:40 and 1:80 configurations. The head-on major merger from Fig.|8]is also presented 
(solid line). 
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Figure 10. Final shape of merger remnant. On the x-axis we show the distance from the core of the merger remnant. The three panels 
present the axis ratios b/a, c/a and b/c (top to bottom, with a > b c). As final measurement location we took 2 X Rh- The results 
for the minor mergers are averaged over all 10 different Plummer realisations. The differences in the cores of the 1:5 and 1:10 mergers is 
caused by the lower resolution with which these models have been run (see text for details). 



Comparing the cold (Q ^ 0.1) and warm (Q > 0.1) con- 
figurations we notice that at the end of the simulation the 
cold configurations have a higher mass than the warm case. 
In the warm configurations the child galaxies have a higher 
velocity which increases their merging times. However, both 
warm and cold show a size increase larger than that of the 
major merger simulations. This is confirmed by the lines 
that show exponential fit^f] through the configurations with 
Q = 0, Q = 0.1 and Q = 1.0. 

The fits indicate that if all child galaxies merge with 
the primary, the merger remnant would have a size that is 
at least a factor of two larger than the remnants formed by 
major merging independent of Q. For the Q — mergers this 
is even three times as large as the major mergers. Also visible 
in Fig.[TT]is the difference between the R = 200 kpc and R = 
500 kpc configurations. Here we notice, as in Fig.JS] that the 
R — 500 kpc configurations show a larger size increase than 
the R = 200 kpc configurations, even though the amount of 
completed mergers in the R = 500 kpc is less than in the 
matching R — 200 kpc configuration. This is attributed by 
the higher energy input that the child galaxies give in the 
R — 500 kpc configuration compared to the R = 200 kpc 
situation where half the child galaxies are already within 
the dark matter halo of the primary galaxy at T = 0. 



5 Note that the fits through the Q = 1.0 results are based on only 
4 (1) merger events in the R = 200 kpc (R = 500 kpc) configura- 
tions. The observed trend however is the same as for the Q = 0.1 
configuration which is based on 10 (9) mergers. 
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Figure 11. The mass, M, of the merger remnant versus the 
half-mass radius, R^, for major merger configurations and a 1:10 
merger configuration with different virial temperatures (Q). The 
minor mergers are all based on the same Plummer distribution, 
but rescaled to a virial radius of R = 200 (R = 500) in the top 
(bottom) panel and a virial temperature which is varied between 
(cold) and 1 (warm). Both panels show the major mergers with 
P = 0, 10 and 50 kpc. The lines are exponential fits through the 
configurations with Q = (solid), Q = 0.1 (striped) and Q = 1.0 
(dotted). With exponents 0.57, 0.49 and 0.485 for R = 200 and 
0.62, 0.49 and 0.484 for R = 500. Configurations are plotted after 
each merger event and therefore occur multiple times. 
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5 DISCUSSION 

The simulations indicate that the effects of the minor merg- 
ers becomes more pronounced when more mergers are in- 
volved in the mass growth. To test if this trend continues 
we performed three extra simulations with 160 child galax- 
ies (mass ratio 1:160). With this many children the simu- 
lation approaches the regime of a continuous stream of in- 
falling material. The three simulations gave very similar re- 
sults with little variation between the runs. With this many 
child galaxies their distribution around the primary galaxy 
becomes almost uniform. The effect of the random place- 
ment of the child galaxies is therefore much smaller than 
in the simulations for 5 to 20 child galaxies. Furthermore 
some of the children will merge with other children before 
they reach the primary galaxy. For the 1:160 simulations we 
use N = 8.8 x 10 6 particles evenly divided over the primary 
(N = 4.4 x 10 6 ) and the child galaxies (N = 27500 per 
galaxy). 



5.1 Properties of the merger remnant 



The circular velocity (V c = y/M(< r)/r ) and the cu- 
mulative mass profile at the end of the simulations (at 
T = lOGyr), both are presented in Fig. 1121 The top panel 
shows V c and the bottom panel presents the cumulative mass 
profile over the same distance from the core. For the mergers 
we averaged the results of the random realizations. 

The major mergers cause an increase in the V c com- 
pared to the isolated model. Whereby the remnants formed 
by major mergers with an elliptic orbit (P = 50kpc) have 
a V c which is ~100km/s higher than that of the isolated 
galaxy. For the remnants that formed by minor mergers we 
see a decrease in V c compared to the isolated galaxy. Espe- 
cially for the 1:40 and 1:80 the effect is quite pronounced, 
with the large size increase comes a decrease in velocity. This 
indicates that the lighter the child galaxies are the lower 
the V c becomes, this trend continuous for the 1:160 mod- 
els. These models have the lowest circular velocity of all the 
configurations while at the same time showing the largest 
size increase. 

The size increase can be deduced from the cumulative 
mass profiles presented in the bottom panel of Fig. 1121 We 
notice that the mass in the core of the merger decreases 
when more child galaxies are involved in the merging pro- 
cess. This effect is especially pronounced for the smaller 
mass-ratios (<1:20). The smaller child galaxies have longer 
merging times, this combined with their interactions inside 
the core of the merger remnant does not allow the merger 
remnant to stabilize. Therefore the mass does not sink back 
to the core within the simulation period. This is also ap- 
parent in the density profiles presented in Fig. 1131 Here the 
major mergers have a density comparable to that of the iso- 
lated galaxy which in turn is higher than the density of the 
minor mergers. For the 1:5 and 1:10 mergers the density 
profiles have a similar shape to that of the isolated galaxy 
and the major merger remnants, but the configurations with 
mass ratios of ^ 1 : 20 have a different shape. The relatively 
sudden depression in the density profiles of 1:80 and 1:160 
near r ~ 1 kpc is a trend that is visible in the 1:40 con- 
figuration and in slightly less pronounced form for the 1:20 
configurations. The longer merging times and dynamic inter- 
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Figure 12. Circular velocity (top) and cumulative mass (bottom) 
profiles after 10 Gyr for different merger configurations. Each line 
shows the averaged result of all random realisations for each of 
the configurations. The velocity and mass profile of the isolated 
model arc indicated with the thick solid line. 
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Figure 13. Density profiles of the different major and minor 
configurations at the end of the simulation (T = 10 Gyr). Also 
displayed is the initial density profile of the isolated galaxy. The 
vertical line is the softening used in the simulations. 



actions that takes place when the number of child galaxies 
is ^ 20 causes a slowdown in circular velocity, a reduction 
in mass and therefore a drop in density in the inner parts of 
the merger remnant. 

The effects described above influence where the mass 
of the child galaxies is accreted in the merger remnant. In 
Fig. [14] we present the fraction of baryonic mass that orig- 
inates from the child galaxies as a function of the distance 
to the center of the merger remnant. The curves are con- 
structed by binning the particles in 1 kpc bins and compute 
the fraction of mass that was brought in by the child galax- 
ies. The curves representing the major mergers (solid-red for 
model P0 and dashed-green for P50) are horizontal and at 
0.5 in Fig. 1141 which indicates that they mix homogeneously. 
Minor mergers experience a different mixing; the precise ef- 
fect depends sensitively on the number of infalling galaxies. 

We take a closer look at the 1:5 to 1:20 mass ratios 
merger s. These configurations behave as expected lHilz et al.l 
l|2012bh : the fraction of child material is small in the core of 
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the merger remnants and high in the outskirts (r ^ 10 kpc). 
While the mass ratio continues to decrease the fraction of 
foreign material in the outskirts also decreases; until the 
majority of outskirt material is original for the most ex- 
treme mass ratio (1:160). We have to note here that for the 
most extreme mass ratio's the merger process has not com- 
pleted at T = lOGyr, which is noticeable in Fig. [14] by the 
mean line appearing smaller than 0.5. In the central portion 
(r 7 kpc) of the most extreme mass ratio mergers the op- 
posite happens, in the sense that the contribution of foreign 
material is actually increasing (in particular noticeable for 
the 1:160 mass ratios). 

In Fig. [15] we present the number density of the for- 
eign and native particles as a function of the radius. We 
do this for a mass ratio of 1:1, 1:10 and 1:160. As we al- 
ready discussed in relation to Fig. [I4]the mass in the major 
merger remnant is evenly distributed between the primary 
and the infalling galaxy. In the 1:10 merger (top panel) the 
two curves for the foreign material and for the native ma- 
terial are clearly separated, in contrast to the lines for the 
major merger. This indicates that the native material is dis- 
tributed differently than the foreign material. The same ef- 
fect is noticeable in the 1:160 merger (lower panel in Fig. ll5[) . 
but less pronounced. This indicates that in the 1:160 merger 
the material of the primary and child galaxies is mixed more 
evenly and approaches a mass distribution similar to that of 
the major mergers, whereas in the 1:10 merger the native 
and foreign material are distributed quite differently. 

The results presented in the previous § indicate that for 
mass ratios 1:1 to 1:20 the accretion behaves according as if 
they grew from inside-out; the material of the child galax- 
ies is st ripped and depos ited in the outside of the primary 
galaxy (|Hilz et al.ll2012bl ). When the number of child galax- 
ies increases beyond 20 this behavior changes in that the 
mixing becomes more homogeneous. This is caused by the 
self-interactions of the child galaxies with the nucleus of the 
primary. These self-interactions prevent the merger galaxy 
from settling to equilibrium and cause the core to heat up. 

We theorize that the continuous bombardment of mi- 
nor galaxies causes the merger remnant to remain quite 
dynamic which reduces the dynamical friction. This in it- 
self allows new incomers to penetrate deeper into the par- 
ent g alaxy. This effect m imics the process of violent relax- 
ation |Lvnden-B cll 1967). This effect becomes more efficient 
when the number of child galaxies increases. During the in- 
fall they interact with other infalling children before they 
break-up, after which their mass is distributed evenly in the 
merger remnant. The interactions are not strong enough to 
permanently eject a child. When the child galaxies are more 
massive 1:5 to 1:20 mergers self interactions between the 
children is negligible and they tend to be disrupted upon 
the first pericenter passage of the merger product. As a re- 
sult their material tends to be deposited in the outskirts. 

The net growth in size observed in our simul ations ex- 
ceeds those report ed in (iNipoti et al.ll2009al . [20121 ). The ini- 
tial conditions of l|Nipoti et al.ll2009a|j2012T l. however, devi- 
ates from ours in the sense that they adopt a more realistic 
mass and shape distribution than in our more theoretical 
approach. As a consequence the results cannot be compared 
trivially. The more efficient growth in our simulations is a 
result in the distribution of the minor galaxies, and not by 
the choice of their densities; the effect of the child densities is 
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Figure 14. Origin of the mass in the merger remnant at T = 
10 Gyr. Along the y-axis we present the fraction of mass that 
originates from the child galaxies at distance r from the center 
of the merger remnant. We binned the mass in 1 kpc bins. When 
the fraction drops to then there is no contribution in mass from 
the child galaxies. 



negligibl e (see App e ndixlBl . Our results ar e consi stent with 
those of iHilz et all [j2012~af ): lOogi fc Habl l|2012h for mass 
ratios < 1 : 20, in the sense that the growth in size is suffi- 
ciently efficient to explain the observational results when the 
number of infalling galaxies > 5 (and consequently a mass 
ratio < 1 : 5). For more extreme mass ratios ^ 1 : 20 we 
find a gradu al change in the be havi or. This regime was n ot 
simulated bv lffilz et all (|2012aT ) and lOoei fc Habd (|2012r i. 



6 CONCLUSION 

We have studied dry minor mergers of galaxies in order to 
understand the observed growth in size without much in- 
creasing the mass of compact massive galaxies. The simula- 
tions ware performed the Bonsai GPU enabled tree code for 
N op to 17.6 million particles. Our simulations start with a 
major galaxy and a number of minor galaxies in a kinemat- 
ically cold (and warm) environment. The total mass of the 
child galaxies equals that of the primary. The simulation 
were performed for 10 Gyr, and we study the size growth 
of the merger remnant. We demonstrate that mergers with 
a mass ratio 1:5 to 1:20 satisfactory explain the observed 
growth in size of compact galaxies. The growth of the merger 
remnant is always at least a factor of two higher than in the 
case of a major merger. This result is robust against varia- 
tions in the initial density of the child galaxies. The mass of 
the minor galaxies tend to be accreted to the outside of the 
merger remnant. As a consequence the core of the merger 
remnant, formed by the original primary galaxy, will hardly 
be affected by the merging process. This 'inside-out-g rowth' 
is consistent with previous studies (|Hilz et al.ll2012bf i. 

If the number of child galaxies exceeds 20 the behav- 
ior of accretion changes, in the sense that the minor galax- 
ies tend to accumulate in the central portion of the merger 
remnant. Evidence for this 'outside-in' growth is present in 
the density and velocity profiles of the merger remnant. The 
outside-in growth is mediated by self interactions among the 
child galaxies before they dissolve in the merger remnant. 
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Figure 15. Particle distribution as function of distance from center, normalized to the volume size. The x-axis shows the distance to 
the remnants core, the y-axis the normalized particle count. This plot is similar as a density distribution, but now split into two lines 
indicating the origin of the particles (mass). The top and bottom panels both show the same averaged major merger result (P0). In 
addition the top panel shows the average result for the 1:10 R = 500 merger and the bottom panel the averaged result for the 1:160 
mergers. Visible is that the 1:10 and major merger show different behaviour, where the major merger shows perfect mixing, the 1:10 
shows a clear distinction between the primary and child galaxy mass distribution. The 1:160 however, shows a distribution that is closer 
to the major merger distribution than to the 1:10 distribution, indicating that 1:160 shows more signs of particle mixing. 



We conclude that the observed large massive elliptical 
galaxies can be evolved from compact galaxies at z £ 2 if 
they have grown in mass by accreting 5-10 minor galaxies 
with a mass ratio of 1:5 to 1:10. The majority of accreted 
mass will be deposited in the outskirts of the merger rem- 
nant. 
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APPENDIX A: RESOLUTION EFFECTS 

As discussed in Sect. [2] we perform a range of simula- 
tions with various resolutions for different major and minor 
merger setups. Depending on the number of child galaxies we 
use low (LR) and high (HR) resolution configurations (see 
Tab.[T}. In the LR configurations we adapted TV = 4.4 x 10 5 
particles of which half belongs to the primary and the other 
half to the child galaxies. For the HR simulations we have 
TV = 4.4 x 10 6 . To verify that the chosen resolutions are 
sufficient we ran a subset of our configurations with even 
higher resolutions. The results can be seen in Fig. IA1I Here 
we present the size evolution of two major merger configu- 
rations, a 1:10 merger with R = 200 kpc and R — 500 kpc 
and one of the 1:80 mergers. Each configuration is run using 
the default resolution, as specified in Tab. [2] and using two 
higher resolutions For the 1:80 simulation this translates to 
TV = 4.4 x 10 6 , TV = 8.8 x 10 6 and TV = 17.2 x 10 6 . For the 
other configurations we use TV = 0.44 x 10 6 , TV = 4.4 x 10 6 
and TV = 8.8 x 10 6 . For each of the resolutions we checked 
if the chosen time-step and softening was sufficient to keep 
the model stable in isolation using the method described 
in Sect. [2] In all but the 1:80 configurations there is little 
to no difference in the size of the merger remnant over the 
course of the simulation. Indicating that our choice of TV is 
sufficiently large. The large amount of child galaxies in the 
1:80 configurations cause child galaxies to interact with each 
other and not directly merge with the primary galaxy. The 
galaxies keep flying in and out of the remnant and makes it 
harder to determine Rh of the merger remnant which causes 
the large fluctuation in size of the 1:80 merger models. How- 
ever, the trend of the Rh increase is independent of the used 
resolution. As final test we performed the 1:80 configuration 
with TV = 4.4 x 10 5 particles (not plotted) which quantita- 
tively gives the same result, indicating that TV = 4.4 x 10 6 
is adequate for the 1:80 configurations. 
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Figure Al. Rh versus time for a subset of configurations. Each 
configuration is run with three different resolutions to test de- 
pendence on the number of particles used. We present the results 
for two major merger configurations (P = and 50kpc), for two 
1:10 minor merger configurations (R = 200 and (R = 500) and 
one 1:80 merger configuration. 



Model 


Bulge mass 


Cut-off 


Rh 


P 




[M ] 


[kpc] 


[kpc] 


[lO lo M /fepc 3 ] 


Primary 


2 x 10 11 


100 


1.36 


1232 


standard 


2 x 10 10 


10 


1.26 


142 


compact 


2 x 10 10 


5.1 


0.68 


307 


supercompact 


2 x 10 10 


0.9 


0.25 


4609 



Table Bl. Properties of the extra 1:10 child configurations. The 
first column indicates the model, either the primary or one of 
the three child configurations. The second column indicates the 
bulge mass. The third the cut-off radius for the dark matter halo, 
which is used to configure the galaxies. The fourth (fifth) column 
indicates the half-mass radius (density) of the galaxy in isolation 
at T=l Gyr. 



APPENDIX B: 
DENSITY 



THE EFFECT OF CHILD 



As described in Sect. [2] we do not change the density of the 
child galaxies when generating the initial conditions. Instead 
we change the cut-off radius to be a factor, /, smaller than 
that of the primary galaxy. This / scales with the mass ratio 
of the child galaxy. For example the cut-off radius of a child 
galaxy with mass ratio 1:10 is 10 times smaller than that 
of the primary. Together with the cut-off radius the mass 
is also reduced by a factor 10 giving the child galaxies a 
density lower than that of the primary galaxy (§01 • 

To test the effect of the density on the size growth of the 
merger product we generated two new 1:10 child galaxies. In 
addition to the standard child galaxy we created a compact 
and a supercompact child galaxy. These galaxies are gen- 
erated such that their densities are higher than that of the 
standard child galaxy. The galaxy properties are presented 
in Tab. [Bl] 

With the new child galaxies in place we selected one of 
the 1:10 R — 500 kpc merger configurations and ran this re- 
alization with the new child galaxies. Each model is evolved 
for 10 Gyr using high resolution (N = 4.4 x 10 6 ). 

The effect of the child density on the size growth of the 
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Figure Bl. Size evolution. The configurations are based on the 
same primary galaxy and Plummer distribution, but the child 
galaxies have different size and density properties. 



merger remnant is presented in Fig. IB1I where we show Rh 
as a function of time. The standard and compact config- 
urations show similar evolution. The supercompact config- 
uration shows a different behaviour during the first 5 Gyr 
of the simulation. This is an artifact of the way we deter- 
mine the size of the merger remnant which does not work 
perfect in this particular configuration (see §[4|) . The com- 
pactness of the galaxies causes the cumulative mass-profile 
to be irregular compared to the simulations with less dense 
child galaxies. As a consequence we are not always able to 
correctly measure the position in the profile that indicates 
the end of the merger remnant. After 5 Gyr the procedure 
works correct, because at this point the child galaxies have 
been involved in so many interactions that the mass-profile 
is smooth again. When the method starts working again we 
notice that all three simulations show similar size growth, 
independent of their density. 



Bl Circular velocity 

The circular velocity of the standard and high density child 
simulations is presented in Fig. IB2I (top) we notice that 
the density of the child galaxy does affect the results. The 
simulation with the most dense child galaxies (striped line 
with stars) has a circular velocity similar to that of the 
isolated galaxy (thick solid line). The compact configura- 
tion has a circular velocity that is ~ 100 km/s lower than 
the supercompact configuration. The standard configura- 
tion has a circular velocity in between the results of the 
compact and supercompact models. There is no clear trend 
in the effect the child density has on the circular velocity. A 
similar effect can be seen in the cumulative mass distribu- 
tion of Fig. IB2I (bottom). The difference between the runs 
is caused by the changed merger history. Even though the 
configurations are based on identical Plummer realisations 
the difference in density causes the merger history to be al- 
tered. Tidal stripping has less effect on the more compact 
galaxies when they move through the dark-matter halo of 
the primary than it has on the standard galaxy. Instead the 
higher density causes these galaxies to stay tightly bound re- 
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Figure B2. Circular velocity (top) and cumulative mass (bot- 
tom) after 10 Gyr for different merger configurations. The merg- 
ers are based on the same Plummer distribution but the child 
galaxies have different densities. The velocity and mass profile of 
the isolated model are indicated with the thick solid line. 
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Figure B3. Origin of the mass in the merger remnant. Shown on 
the y-axis the fraction of mass that originates from the primary 
galaxy at a certain radius from the center. The mass is binned in 
1 kpc sized bins. If the fraction is then there is no contribution 
in mass from the child galaxies. 



suiting in interactions between other children and the core 
of the primary galaxy. Because of these interactions some of 
the child galaxies will be flung from the galaxy during the 
first passage instead of merging with the primary. 

In Fig. IB3l we present the distribution of accreted mate- 
rial throughout the merger remnant. In this figure we notice 
large differences between the different configurations caused 
by the differences in child density. The solid line represents 
the standard simulation where the mass is accreted on the 
outside of the primary (outside- in formation). For the denser 
children the mass is accreted closer to the core of the primary 
galaxy. This is the result of the less effective tidal-stripping 
described in the previous §. The dense child galaxies are able 
to progress further to the core of the merger remnant before 
they break up and are accreted. While the standard galax- 
ies loose much of their material outside the core. This effect 
is best indicated by the supercompact configuration. In the 
merger remnant of this simulation most of the mass in the 
inner 8 kpc originates from the child galaxies. While in the 
other configurations most of the mass in this inner part of 
the merger remnant originates from the primary. 

Concluding we can say that the density has no effect on 
the size of the merger remnant, but it does have an effect 
on how the remnant is build-up. With either mass being 
accreted on the outskirts (low density children) or in the 
inner parts of the merger remnant (high density children). 



